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This talk discusses some recent studies of gluon saturation in nuclei. We stress the connec- 
tion between the initial condition in heavy ion collisions and observables in deep inelastic 
scattering (DIS). The dominant degree of freedom in the small x nuclear wavefunction is 
a nonperturbatively strong classical gluon field, which determines the initial condition for 
the glasma fields in the initial stages of a heavy ion collision. A correlator of Wilson lines 
from the same classical fields, known as the dipole cross section, can be used to compute 
many inclusive and exclusive observables in DIS. 



1 Connection between small x DIS and HIC: Wilson line 

The initial condition in a heavy ion collision (HIC) is determined by the wavefunctions of 
the two colliding nuclei, parametrized by and x. As in any hadronic collision, the typical 
magnitudes of these parameters can be estimated as ^ (P-l)'^ ^-iid x ^ {p±)/^, where {p±) 
is the typical transverse momentum of the particles being produced, and ^/s is the collision 
energy. At rclativistic energies, such as at RHIC and LHC, this means that the relevant domain 
for bulk particle production is at very small x. Gluon brehmsstrahlung processes lead to an 
exponentially (in rapidity y — Inl/x) growing cascade of gluons in the wavefunction. The 
number of gluons in the wavefunction grows as dN/ dy ~ a;~^, where the phenomenologically 
observed value is A ^ 0.2 . . .0.3. When the number of gluons grows large enough, eventually 
their phase space density becomes large, with occupation numbers ^ l/ckg; in terms of the field 
strength this meas ^ 1/g. At this point the nonlinear terms in the QCD Lagrangian (think 
of the two terms in the covariant derivative 9^ + igA^ ) become of the same order as the linear 
ones, and the dynamics becomes nonperturbative. Due to the nonlinear interactions the gluon 
number cannot grow indefinitely, but it must saturate at some l/ofs for gluons with p± < Qg, 
where Qg is the saturation scale. When ^ x~^ becomes large enough, as{Qs) ^ 1 and the 
dynamics of these fields is classical. This situation is most conveniently described using the 
effective theory known as the Color Glass Condensate (CGC, [1]), where the large x degrees of 
freedom are described as a classical color current and the small x gluons as classical fields 
that this current radiates: [Dfj,,F'^'^] — J". 

Let us consider a hadron or a nucleus moving in the +z-direction. Its color current in 
the CGC formalism has only one large component, the one in the +-direction (recall that 
x^ = (t±z)/'\/2) . For a nucleus moving at high energy we can take the current to be independent 
of the light cone time x~ as J"*" = p{xt,x~) with a very narrow, (5-functionlike support in x~: 
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Figure 1: Left: Spacetime structure of the field. In regions (1) and (2) there is a transverse 
pure gauge field Q with no energy density. In region (3) after the collision there is the glasma 
field. Right: 



p{'x.T,x~) ^ (5(x~)p(xt). a simple solution for the equations of motion can be found as 
= p{:x.T, a;~)/V^; this is known as the covariant gauge solution. In order to have a physical 
partonic picture of the gluonic degrees of freedom it is necessary to gauge transform this solution 
to the light cone gauge = 0. The gauge transformation that achieves this is done with the 
Wilson line constructed from this gauge field: 



C/(xt) = Pcxp 



(1) 



This results in a field with only transverse components: 



^(1,2) -f^(l,2)(XT)5.t/(V2)(XT) (2) 

for both of the colliding nuclei (1,2) separately. The intial condition for the "glasma" fields 
at r > is given in terms of these pure gauges [3] . 

Inside the future light cone r > the field equations must be solved either numerically or in 
some approximation scheme. The spacetime structure described here is illustrated in Fig. [T] 
(left). In the rest of this talk we shall be referring to the numerical "CYM" (Classical Yang- 
Mills) computations [4]. 

To see the connection to DIS it is convenient to consider the process in a Lorentz frame 
where the virtual photon has a large longitudinal momentum. In the target rest frame (or more 
properly the "dipole frame" [S] that does not leave all the high energy evolution in the probe) 
the timescales of the quantum fluctuations of the virtual photon are extremely slow. In order to 
interact with a hadronic target it must therefore split into a quark-antiquark pair already long 
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before the scattering. This qq-dipo\e then interacts with the hadronic target with a scattering 
ampHtude whose imagimary part is known as the "dipolc cross section" . As typical hadronic 
scattering ampHtudes at high energy, that of the dipolc is almost purely imaginary, and we shall 
here neglect the real part. The dipole cross section can be obtained from the quark propagator 
in the gluonic background field of the target, which is quite naturally given by the same Wilson 
line 111) 6 . The dipole cross section (which, in general, is a function of the size of the dipole 
Tt, the impact parameter and x) is the correlator of two Wilson lines 



aivT) = I d^bT^ (l - C/t (br + y) (br - Y 



(4) 



For example, the total virtual photon cross section can be obtained by convoluting the dipole 
cross section with the virtual photon wavefunction which relates the of the photon to the 
size of the dipole r ^ l/Q- 



-7 P 

'l,t 



d hr d tt dz "^I riQ ^^t,z) (Tdip(a;, tt, br) 



(5) 



Fourier-transforming instead of simply integrating over the impact parameter dependence gives 
access to the momentum transfer to the target in diffractive scattering. The inclusive diffractive 
virtual photon cross section (really the elastic dipole-photon cross section) is proportional to 
the square of the dipole cross section 



D,tot 
^L.T 

dt 



1 

16^ 



d^TT 



^lj,{Q'^,rT,z) a^ip(a;,rT, A) 



(6) 



and diffractive vector meson production can be obtained by projecting on the virtual photon 
wavefunction 



D,V 

dt 



1 

16^ 



d^TT / dz U-^^- 



L,T 



((5^,rT,2:)crdip(a;,rT, A) 



(7) 



The exclusive cross sections are proportional to the dipole cross section (the scattering am- 
plitude) squared, whereas the inclusive one depends on it linearly; this is due to the optical 
theorem and our approximation that the scattering amplitude is purely imaginary. We shall 
now go on to discuss some recent applications of saturation ideas to heavy ion collisions and 
DIS phenomenology, trying to stress the unity of the approach between the two. 



2 Gluon multiplicity at RHIC and LHC AA collisions 

Ideally one would like to measure the value of Qs in DIS experiments and use the resulting 
value as an independent input in calculations of the initial state of heavy ion collisions. In 
practice most of the exiting CYM computations of the giasma fields have been performed in 
the MV model [9] in terms of the color charge density parameter g^^ that parametrizes the 
fluctuations of the classical color currents J^. One must therefore relate the values of g^/i and 
Qs in a consistent way. In practivc this can be done by computing the Wilson line correlator in 
the MV model, using exactly the same numerical implementation of the model as in the CYM 
calculations, and extracting the correlation length 1/Qs [10] . 
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Figure 2: Left: Comparison of the fit [7] to existing nuclear DIS data from NMC. Right: 
Extrapolation of the gluon multiplicity to LHC energies, from [8] 



The other ingredient necessary in using the existing DIS data to calculate initial conditions 
for heavy ion collisions is the correct implementation of the nuclear geometry in extending 
the paramctrization from protons to nuclei. In a "Glauber" -like formulation of essentially 
independent scatterings of the dipole on each of the nuclcons this is a straightforward estimate, 
see e.g. Refs. [71[TT]. A simple geometrical argument would give the estimate « O.bQ'^pA^^^ , 
where the coefficient in front follows from the internucleon distance in a nucleus being smaller 
than the nucleon radius. The actual values in the estimate of Ref. [7] are shown in Fig. [1] (right). 
For other estimates of Qs based on DIS data see Ref. [12]. Being really consistent with high 
energy evolution would require some further theoretical advances, since the approximation of 
independent dipole-nucleon scatterings will break down during the evolution. In the infinite 
momentum frame this can be thought of as gluons from different nuclcons starting to interact 
with each other. 

Combining these ingredients the CYM calculations [4] of gluon production paint a fairly 
consistent picture of gluon production at RHIC energies. The estimated value Qs ~ 1-2 GeV 
from HERA data [71 [Tl] (corresponding to the MV model parameter g^n sa 2.1 GeV [TUj) gives 
a good description of existing nuclear DIS data from the NMC collaboration, see Fig. [2] (left). 
The same value leads to « 1100 gluons in the initial stage. Assuming a rapid thermalization 
and nearly ideal hydrodynamical evolution this is consistent with the observed ~ 700 charged 
(~ 1100 total) particles produced in a unit of rapidity in central collisions. 

The gluon multiplicity is, across different parametrizations, to a very good approximation 
proportional to TiI{\Q^/as. Thus the predictions for LHC collisions depend mostly on the 
energy dependence of Qs- On this front there is perhaps more uncertainty than is generally 
acknowledged, the estimates for A = dlnQg/dlnl/x varying between A = 0.29 [13] and A = 
0.18 [13] in fixed coupling fits to HERA data, with a running coupling solution of the BK 
equation giving something in between these values }15j. This dominates the uncertainty in 
predictions for the LHC multiplicity, see Fig. [2] (right). 
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3 Multiplicity distributions 



One very recent application of the CGC framework has been computing the probabihty distri- 
bution of the number of gluons in the glasma |16| . The dominant contributions to multiparticle 
correlations come from diagrams that are disconnected for fixed sources and become connected 
only after averaging over the color charge configurations. In other words, the dominant correla- 
tions are those arising from resummcd large logarithms of the collision energy and are present 
already in the initial wavefunctions of the colliding nuclei. 

Working with the MV model Gaussian probability distribution 



W[p] = exp 



^2^^ P''(xt)p"(xt) 



(8) 



computing the correlations in the linearized approximation is a simple combinatorial problem. 
The result can be expressed in terms of two parameters, the mean multiplicity n, and a param- 
eter k describing the width of the distribution. The g'th factorial moment of the multiplicity 
ditribution is, to leading order in as, proportional to 2''{q — 1)!. Explicitly, the connected parts 
of the moments = {N'^) are 



(q-l)!fc(^-j with (9) 



27r 



(10) 



n = fN—QlS^. (11) 



as 

These moments define a negative binomial distribution with parameters k and n, which has been 
used as a phenomenological observation in high energy hadron and nuclear collisions already for 
a long time |17| . In terms of the glasma flux tube picture this result has a natural interpretation. 
The transverse area of a typical flux tube is 1/Ql, and thus there are Q^S± ~ iVpT independent 
ones. Each of these radiates particles independently into N^^ — 1 color states in a Bose-Einstein 
distribution (see e.g. [H]). A sum of A: w Nyt{Nc^ — 1) independent Bose-Einstein-distributions 
is precisely equivalent to a negative binomial distribution with parameter k. 



4 Inclusive nuclear diffraction at eRHIC and LHeC 

The large fraction of diffractive events observed at HERA shows that modern colliders are 
approaching the nonlinear regime of QCD, where gluon saturation and unitarization effects 
become important. It should be possible to perform the same measurements in DIS off nuclei. 
There are plans for several facilities capable of high energy nuclear DIS experiments, as the 
EIC [T^ and LHeC [20] colliders. Due to the difficulty in measuring an intact recoil nucleus 
deflected by a small angle, diffractive eA collisions present an experimental challenge. But if they 
are successful, nuclear diffractive DIS (DDIS) would provide a good test of our understanding of 
high energy QCD. Measuring the momentum transfer t in both coherent (nucleus stays intact) 
and incoherent (nucleus breaks up into nucleons) would enable one to go measure directly 
the transverse structure of the gluonic degrees of freedom [5T] instead of the electric charge 
distribution that is measured in low energy experiments. Figure [3] (left) demonstrates some 
expected results from such a measurement. The diffractive structure function can be divided 
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Figure 3: Left: nuclear modification of the diffractive structure function, showing characteristic 
suppression at small /3 (large mass of the diffractive system) and enhancement at large p. From 
Ref. [22| . Right: dominant impact parameters for the different contributions to the proton 
diffractive structure function. 



into different components according to the polarization state of the virtual photon and the 
inclusion of higner Fock states (e.g. qqg in addition to qq) in the dipole wavefunction. All of 
these have different dependences on the impact parameter of the dipole-target collision (see 
Fig. [3] right), which stresses the importance of having a detailed picture of the transverse 
geometry of both the proton and the nucleus. 
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